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Abstract 



We study a two parameter single trace 3-matrix model with SO(3) global symmetry. The model has two 
phases, a fuzzy sphere phase and a matrix phase. Configurations in the matrix phase are consistent with 
fluctuations around a background of commuting matrices whose eigenvalues are confined to the interior of 
a ball of radius R = 2.0. We study the co-existence curve of the model and find evidence that it has two 
distinct portions one with a discontinuous internal energy yet critical fluctuations of the specific heat but 
only on the low temperature side of the transition and the other portion has a continuous internal energy 
with a discontinuous specific heat of finite jump. We study in detail the eigenvalue distributions of different 
,1^ ' observables. 

1 Introduction 

One of the most striking and interesting features of multi-matrix models is the phenomenon of emergent 
geometry. The notion of classical geometry changes drastically within the context of matrix models; the 
geometry is no longer a basic concept which exists a priori but instead it emerges dynamically as a conse- 
quence of the reordering of degrees of freedom. This is in many ways exciting because matrix models can 
lead to new ways of thinking about the structure of the space time. 

The interest in matrix models has grown since it was suggested that they might provide a non-perturbative 
definition for M theory [T]. Several kinds of matrix models have been proposed for this purpose [TJ[2j[3]. 
The IIB matrix model (IKKT model) [2] is one of these proposals; it is a large N reduced model [4] of 
ten-dimensional supersymmetric Yang-Mills theory and the action is a matrix regularised form of the Green- 
Schwarz action of the IIB superstring. It is postulated that it gives a constructive definition of type IIB 
superstring theory. 

Finite dimensional matrix models have also been used to regulate field theories [7J [8] [9] 1101 1111 1121 113] 
and diverse fully non-perturbative numerical studies have been performed, see for instance |14l 1151 1161 1171 

□sunn]- 

The pure commutator action, which we refer to as the Yang-Mills matrix model, is also particularly 
interesting since in d = 10, 6 and 4 dimensions it corresponds to the bosonic part of the IKKT model. The 
model is well defined in dimensions d > 2 and matrix size N > 3 [29]. In dimension d — 2 the model, with a 
quadratic term is added to stabilise it, is exactly solvable 24, 32]. These bosonic Yang-Mills matrix models 
in different dimensions have been considered as possible realizations of emergent geometry/gravity [30] . 

Numerical studies of pure Yang-Mills matrix models were performed in [4] for different dimensions (num- 
bers of matrices) and more recently in [5j [6] it is argued that perturbation theory around a background 
of commuting matrices gives a good approximation to the 3-dimensional model. In [6] based on a two 
loop computation it is predicted that the eigenvalues of the background commuting matrices are uniformly 
distributed within a ball of radius R ~ 1.8 which is in broad agreement with our findings here. 

Our starting point is an action in which the basic objects are simple Hermitian matrices; no geometrical 
background is assumed a priori. The model describes the statistical fluctuations of matrices with prescribed 
energy functional. The geometry arises as a condensate around which the system fluctuates. Generically, 
multi-matrix models can undergo transitions between different geometries and phases with no geometrical 
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content. Given the novelty of these phenomena it is worth studying the simplest model that exhibits such 
phenomena in detail. 

The simplest model in which a geometry has been shown to be emergent is a 3-matrix model consisting 
of the trace of the square of the commutator of the matrices (a Yang-Mills term) plus the epsilon-tensor 
contracted with the trace of the three matrices (a Myers term) [22] . This model was studied in |15lll7|[l9ll20| 
and is a static bosonic subsector of the BMN model [3]. 

The model exhibits a geometrical phase for sufficiently large coupling to the cubic Myers term, with the 
geometry being that of a fuzzy sphere [231 124] ; a non-commutative [!TT] version of the commutative sphere 

02]. 

At a critical coupling, which can be traded for a critical temperature, a phase transition occurs and the 
condensed geometry evaporates. In the geometrical, low temperature phase, small fluctuations around this 
condensate correspond to a U(l) gauge and scalar field multiplet 13, 251 I2U]. 

The model can be extended by adding appropriate potentials in order to enhance the range of parameters 
in which the fuzzy sphere phase is stable [20] . The general phase diagram for a class of such models was 
predicted in [20] , one of the purposes of this paper is to check these predictions in detain in a non-perturbative 
study. We find that indeed the phase diagram is well predicted by theoretical expressions presented in [20] . 

Manifolds in higher dimensions can emerge from more general matrix models with essentially the same 
structure and phenomenology as can be seen in [261 1271 [28] . 

In the current paper we study, in detail, the high temperature phase of the model and establish that this 
phase is consistent with the three matrices fluctuating around commuting matrices where the eigenvalues of 
the commuting matrices are confined to the interior of a solid three dimensional ball. The fluctuations are 
still significant as the commutator of pairs of matrices is itself a peaked (almost triangular) distribution. 

We concentrate on the effect of adding a quadratic, mass like perturbation to the 3-matrix model. We find 
that below a critical temperature any negative massive perturbation induces a transition from the matrix 
phase to the fuzzy sphere phase. For higher temperatures a more negative quadratic coupling is necessary. 
The phase diagram, which is one of the principal results of this paper, is shown in figure 3. 

• We find the ground states of the system are characterised by either the N dimensional irreducible 
representation of 5(7(2) (the fuzzy sphere phase) or a continuum spectrum (matrix phase) for one of 
the matrices. These characterize the two phases of the system. Though meta stable states other than 
the irreducible fuzzy sphere are present in the system they never correspond to the true ground state, 
they do however appear stable in the large N limit where tunneling is suppressed. 

• We study the phase diagram (see figure 3) as a function of the two parameters r and a and locate the 
co-existence curve with some precision. 

• The co-existence curve rapidly asymptotes to the special line r = § where the energy functional defining 
the model becomes a complete square. 

• We found evidence for two distinct types of transition in the emergent geometry scenario: for < r < | , 
as the transition is approached from the fuzzy sphere phase with fixed r, the system has a divergent 
specific heat with critical exponent a = | while crossing the phase boundary for fixed a > 4.02 there 
appear to be no critical fluctuations and the transition is one with a continuous internal energy and 
discontinuous specific heat. 

• We find that a useful description of the matrix phase is in terms of fluctuations about a background of 
commuting matrices whose eigenvalues are concentrated in a three dimensional ball of radius R — 2.0. 

The paper is structured as follows: Section 2 describes the model and the predictions of the phase diagram 
from the effective potential [20] in the fuzzy sphere phase. Section 3 describes the matrix phase and the 
consequences of a background of commuting matrices with eigenvalues uniformly distributed within a ball. 
Section 4 describes our numerical results and section 5 contains our conclusions. 

2 The model 

We begin by considering an action functional built as a single trace quartic polynomial of 3 Hermitian 
matrices with SO(3) symmetry. There are four available invariants: 

Tr D a D b D a D b , Tr {D 2 a f , Tr ie abc D a D b D c , Tr D\. (1) 



In this paper we restrict our study to the two parameter model given by the action 

1, 



S[D] = f Tr 



-j [D a , D b ] 2 + - e abc [D a D b ] D c + TD 2 a 



(2) 



where stability of the model requires the r > 0. For r < the action for commuting matrices is unbounded 
from below. The most general model would include in addition Tr(D 2 ) 2 with an extra coupling. 

The action ([5} is invariant under unitary transformations U(N), D a — ¥ UD a U^ and global 5*0(3) rota- 
tions of the matrices. The modes c a = Tr(D a ) decouple from the others and we therefore choose to work 
with traceless matrices, Tr D a = 0. The parameters of the model are r and a. a can be identified either as 
the Yang-Mills coupling constant g 2 = a~ 4 or the temperature T = a -4 [20] , 

Saddle points of the action, derived from the condition 8S — 0, are given by solutions of 

[D b ,iF ab ]+2rD a =0, with F ab = i[D a , D b ] + e abc D c 

and include the trivial solution D a = and D a = 4>Ja where J a are representations of SU(2), not-necessarily 
irreducible. For D a = (j>J a , then F ab = — (<f) 2 — 4>)e abc J c we get an algebraic equation for 4> given by 



:3 ,2 



+ t<j>) = (3) 



t = {^-,<P + } = <0, 1 -±^E, 1 —4EE\. (4) 



the explicit solutions are 



The first corresponds to D a — and is the ground state of the system for r > i , is a local minimum for 
2 > t > and though it appears to be a local maximum for r < 0, the model has no ground state for such 
value of r. The second solution gives the local minima D a = <fr-J a , for r < \. These minima D a — and 
D a — 4>-J a are separated by a potential barrier whose highest point is localised at the local maximum </>+. 
As r approaches zero <j>- —¥ 1, which corresponds to the case for rn 2 — r = studied in |15l 1191 120| . For 
arbitrary r has been also studied by other authors [211120] , 

For the configurations D a — 4>Ja the action 

S[D] = V clas M) = ^][>a,(m) & A (^ - y + ^ ) (5) 

where Ci(ni) is the Casimir of the representation of dimension m and J^i n » = N. For r < | the potential 
14; a ss(0-) < and the configurations D a = 4>-Ja is of lower energy than D a — 0, which is zero. The 
configuration with minimum energy is then given by maximising the sum of Casimirs, this is achieved by 
J a — L a , where L a is the irreducible representation of dimension N. For r = | we have V c iass (0-) = and 
the configurations D a = and D a — cf>-J a , with any representation J a , become degenerate. Figure [2] shows 
the potential ([5]) for different r. 

Therefore the classical prediction is that for r < |, and any value of a the ground state is D a — 4>-L a 
and small fluctuations around this configuration have the geometrical content of a Yang-Mills and scalar 
multiplet on a background fuzzy sphere. 

The parameter domain | < r < \ is of special interest. The classical analysis above suggests that the 
fuzzy sphere is unstable here, however, since the potential still has a local minimum and in the TV — >■ oo 
limit all representations such that -k £V mCiini) —¥ become degenerate with the D a = configuration 
one might wonder if one of these configurations gives the ground state of the system when fluctuations are 
included or perhaps as suggested by the analysis of 20 and discussed in the next section, the fuzzy sphere 
phase is stabilised by fluctuation. The answer as we will see, is that numerical simulations do not support 
the assertion that the fuzzy sphere is stabilised in this parameter range, but rather that the coexistence 
curve between the fuzzy sphere phase and the matrix phase asymptotes to the line r = |. 

For r = § the model is indeed special, the action is always positive semi-definite and can be written in 
the form 

5 4 „ fi r „ „ , 1 



S[D] = —Tri-[D a ,D b ] + -e abc D c ) , (6) 



from which we see that D a — and the orbit D a = U%L a W have zero action 
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2.1 Quantum corrections and critical behaviour 

In the previous section we described the classical potential predictions for the phase transition by choosing 
D a = <j>L a and studying the potential for <f>. This analysis suggests that (j> plays the role of order parameter 
for the transitions which are taking place in the model. We will now take into account the quantum 
fluctuations. The computation of the quantum effective action and the quantum effective potential was 
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Figure 1: Classical potential for different values of r for fixed a = 5. 



carried out in detail [251 119] for a more general SO(3) globally invariant matrix model. The results there 
correspond to the present case by setting m 2 = 0. 

Using the standard background field method around the classical configuration D — 4>L a one finds the 
effective potential in the large N limit is given by 



2C 2 



T-y + r y 



+ log< 



(7) 



As discussed in [251 119] the phase diagram including fluctuations is obtained from the minimum of this 
effective potential. Due to the log term the effective potential ([7]) is not bounded from below near (j> — and 
since in its derivation it was assumed that we were expanding around D a — 4>L a it is only valid where such 
a ground state exists. 

The minimum is therefore one of the roots of the polynomial <f> — °I ! } — which gives the equation 

4> 4 - 3 + rtf + A = 0, (8) 



and determines 4> in D a — cf>L a at quantum level. Explicitly the minimum is given by 



1 1/1 2r , 1 

4 + 2Vl~T +d+ 2 
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+ d 



(9) 



with the definitions 
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(10) 
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Figure 2: (a) Effective potential for different values of r for a = 5. (b) Effective potential for r = 0.18 for 
different values of a 



V c s is plotted in figure 2. In figure 2(a) the potential is shown for a fixed value of a — 5 for different 
values of r. In figure 2(b) we fix the value of r = 0.18 varying a. From both pictures we can see that there 
exists a region in which the fuzzy sphere ceases to exist. This means that there exists a critical coexistence 
curve in the (r, a) plane where the model undergoes a phase transition between the two phases. 

These expressions predict that for sufficiently large a or low temperature and r < -j the fuzzy sphere is 
the ground state of the model. The local minimum disappears for the simultaneous solution of V^g((j>) = 
and KffO) = so that 



Acj)i - 3(j>, + 2t = 



(11) 



which gives the critical values for <fi and a. 



32r 



and 



2r) 



s 



(12) 



These equations defines the phase diagram (r, a) for the present model. They prediction that at zero 
temperature i.e. a = oo where (j>* — 1/2 and r = 1/4. 

The phase diagram is presented in Figure 3 together with a comparison with numerical simulations. 
Numerical results are in excellent agreement with the prediction from effective potential for r < § where 



the coexistence curve predicted in (|12[1 clearly delimits the fuzzy sphere phase from the matrix phase. The 
fuzzy sphere exists for a> a* and for t < t*, where T,(a) is obtained by inverting (|12l) . 



However, for 



> t > 



simulations show that the coexistence curve asymptotes rapidly to the line r : 



/ \ 1/4 

for a greater than the special value a* = 12 I * r — ) ~ 4.02. As we cross the critical value of a, a 

rather exotic phase transition occurs where the geometry disappears as the temperature is increased. In 
the high temperature phase, which we call a matrix phase, the order parameter (f> goes to zero as N~ . In 
this phase the fluctuations are insensitive to the value of a, they are in fact fluctuations around commuting 
matrices. Since, in the large iV limit, a is unimportant we can rescale the matrices to eliminate a from the 
quadratic term, defining X a — -j=D a , we obtain 



S[X] = NTr(-hx a ,X b ] 2 + ^=e abc X a X b X c + °~ T 
4 3v N 



N 



xl) 



(13) 



and see that in this rescaled model both r and a drop out of the model in the large N limit. 
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Figure 3: Phase diagram (r, a). The solid line corresponds to the theoretical prediction eq. (fl"2j) for r < r* = §. 
The numerical points are obtained from simulations with matrix sizes N — 19, N — 24, N = 35 and TV = 64. 
We see some rounding of the coexistence curve near (a* = 4.02, r = |), before it asymptotes to the line r = |. 



Next we discuss the numerical results which include the eigenvalue distributions for the configurations. 



2.2 Specific heat and phase transitions. 

There are several useful general identities that we can derive between expectation values of observables, the 
simplest of these is the relation between different components of the action. Since our action is a polynomial 
in the matrices X a it can be expressed as S — $4 + S3 + & where Sfc[AX] = X S[X]. Then if we scale the 
fields X a — > XX a in both the action and the measure, the partition function Z, is invariant and under an 
infinitesimal rescaling we obtain the constraint 

4 < Sa > +3 < S3 > +2 < S 2 >= 3(/V 2 - 1), (14) 

which allows us to eliminate < £4 >. Defining Sk = < J& > we have for large N obtain 

5 = J + J 53 + \ S *' ( 15 ) 



For our model, with 



we have 



Xa = -4=J>«, (16) 

V N 



Si = -^Tr([X a ,X b } 2 ), (17) 

S 3 = N^iae a bcTr(X a X b X c ), (18) 

S 2 = Na 2 r Tr(X 2 ). (19) 

Where a = -S= . In the fuzzy sphere phase we know fluctuations are around the configuration D a — 4>L a 
and we then have the predictions (see [20] for details) the partition function Z = e _F with the free energy 
F given by 

F(a,r) = ^+31na. (20) 
where Vis is evaluated with <f> given by ©. From this we can then identify 

< S 2 > rSV < S3 > «V 3 < S 4 > 3 a 



+ ^~ (21) 



N 2 4 ' N 2 6 iV 2 4 

In particular for the average of the action S and the specific heat C v — -jk? < (S— < S >) 2 > we have 

„ 3 a 4 T(j} 2 a 4 (j> 3 , 3 a 5 iri „ ,dd> 

S = - H — and C v = - H 0((i - 2r) — . (22) 

4 8 24 4 32 rvr ; d5 ^ ; 

3 The Matrix Phase 

In terms of X a of (|16[) the model ([2]) becomes 

S[X] = TV" Tr 



1,,, „ ,j . i a r „ ,,,,,. q ,, 2 



[X a ,Xi,] + — — j=e a bc[X a Xb\Xc + r-r=rX a 
4 3 y/ N i\ 



(23) 



We see that at high temperature where a — > 0, provided fluctuations in X a do not grow too rapidly with 
either a or TV then in the large TV limit the model reduces to the pure Yang-Mills term 

S[X]=-^Tr[X a ,X b ] 2 . (24) 

Our numerical results support this and we find that < r N a > and ^ae abc < Tr(X a , X b X c ) > are indepen- 
dent of N, see Figures 4. 

We see immediately from (|15p and the reduction to (|24|) that in the high temperature, matrix phase 
we should expect both S — 4 = C v irrespective of the value of a and r. This is in accordance with our 
numerical results shown in Figures 5 and 7 and the earlier results of J20 , but there is some discrepancy with 
Figure 6 which we believe is due to finite size effects as discussed below. 

It is argued in [5] [6] that fluctuations are about a background of commuting matrices whose eigenvalues 
are uniformly distributed in the interior of a solid ball in R . This in turn predicts that the eigenvalue 
distribution of a single matrix, say X3 , is the parabolic distribution 

p{x) = -^(If - x 2 ). (25) 





<s 3 > 



Figure 4: ^tttt as a function of a and jyl as a function of a for r = 0.23, 0.24. We observe collapsed data in 
both plots in the range showed. 



In [6], based on a two-loop approximation, it was estimated that R ~ 1.8. This is in reasonable agreement 
with our numerical results (see Fig 13) which give R = 2.0. In Figure 13 we can see that the numerical 
results for the eigenvalue distribution of X 3 are very well fit by the parabola (|25l) . 

In the same figure we present the distribution for the commutator i[Xi,X2] showing a symmetric dis- 
tribution with support lying in the interval [—1.5, 1.5]. This indicates, that the fluctuations around this 
background are non-commuting variables. 

If we accept the parabolic distribution (|25p and radius R = 2 we can predict 



Therefore we have 



<£(*2)>= 



5 2 = 



R a SR 2 12 

dxip(xi)x 1 = — — = — . 

R 



2 < Tr Xl > a 2 12 

QT N = r 7v"Y' 



(26) 



(27) 



which goes to zero for fixed a in the large N limit. Our numerical results are in very good agreement with 
this result for S2, as we can see in figure 4 (right). The value of the radius extracted from here is R = 2.01, 
in accordance with the value obtained from the eigenvalues of X3. However, as a increases we can observe 
small deviations. The value R = 2 corresponds to the limit of infinite temperature and large N. 

From our numerical simulations we see that S3 goes to zero even more rapidly than 52- In figure 4 we 
see that the curves for different N collapse when we 

<S 3 > / 2 



ATl/2 



= ( -i&£ a bc Tr XaXbXc 



(28) 



mation becomes poorer as r — ; 
data for r = 0.23 is well fit by 

<S 3 > 



- T -»v- find -# behaves linear in a, however the approxi- 

i and higher corrections subleading in N should be taken into account. The 
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JVl/2 

0.24 we find 

<S 3 > 
jyi/2 



0.142533 a 



0.0243563 a 



<S 2 > 

N 2 



<S 2 > 

N 2 



0.557965 a 



0.574869 ol 
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Figure 5: The specific heat C v as a function of a for different r values, r < 2/9. The curve dashed line 
corresponds to expression (|22| for the specific heat. 
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Figure 6: Observables < S >,G V for t — 0.25. The range of a shown corresponds to a < 3.7 



From which we deduce that S3 — » as N 5 . 

Taking this into account we can write the expectation value of the action *%? a > given by eq. (|15l) 

5 = t + -T7t(0-142533 &- 0.799823 a) + r%-^ 

4 4Af l ' N 5 

= T + -T777T7 (0-142533 a) - i(0.799823 a) + ra%\ for r = 0.23, 

4 4iV 1 / 2 4 5 

5 = T + TT7 (0-024356 6 -0.757264 5 2 ) + r^-- 
4 4N K ' N 5 

= ~ + »!i /■, (0.024356 a) - -(0.757264 a 2 ) + to 2 -; for r = 0.24, 

4 47V 1 /2 4 5 

5 = T + -T777T7 (-0-159707 a) - t(0-701843 a) + tu%; for r = 0.25. 

4 47V 1 /2 '4 5 

The expectation values for the action S are shown in figure 8 for these last two particular values of r. 



(31) 

(32) 
(33) 



4 Numerical Results 

We now describe the results of the numerical results for the expectation value of the action < S > and the 
specific heat C v /N 2 —< (S— < S >) 2 > as function of a for a fixed value of r. 

For all values of r < 2/9 we observed that there is a discontinuity in the expectation value of the action 
< S > across the transition line in Figure 3. This jump in S is well predicted by the theoretical expression 



AS 



-4/ 

-a (- 



+ 



24" 



(34) 



with $* given by equation ([9]) evaluated on the transition. The specific heat is non-analytic along the 
transition, but the nature of its non-analyticity differs on the critical line r = I and the critical curve 



d=5.00 N=24 ^-^- 




J0» 




•« 




.-■' 




,.x" 




«■■'' 


- 


1 

0.8 
0.6 
0.4 
0.2 

-0.2 




■ 


2/9 


- 




0.2 2/9 0.25 



0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 



Figure 7: Observables < S >, C v as function of r for a = 5. The transition occurs at r = 0.219 ± 0.005. 



(r < |, 5*(r)). In the former case the specific heat has the standard jump of a 1st order transition while on 
the latter it diverges as the transition is approached from the low temperature or fuzzy sphere phase side of 
the transition but there is no observable increase in the specific heat as the transition is approached from 
the high temperature side. 

To identify the transition we choose that value of a that gives the specific heat and identify this as a,. 
In figure 5 we show the curve of the specific heat C v for fixed values of r = 0.08 and r = 0.15. The critical 
value of the coupling a* measured from simulation is 2.40 ± 0.05 for r = 0.08 whereas we found a value of 
5* = 2.85 ± 0.05 when r = 0.15. These are to be compared to the prediction (|12l) from which we get the 
values 5* = 2.366 and 5* = 2.790 respectively. As we can see the numerical results are in good agreement 
with the our theoretical predictions. 

Thought the transition is characterised by a latent heat, as it was pointed out in [191I20J . it is unusual 
in that there are also critical fluctuations in the specific heat as the transition is approached from the low 
temperature phase for r < |. These lead to a divergent specific heat with specific heat exponent a — 1/2. 
Our simulations are consistent with the expected divergence in the low temperature phase. 

In figure 6, where we show 5 for r = 1/4 and a < 3.7. We see the data collapses when plotted against 
a. We also show that the specific heat is fluctuating around C„ = 4. These are consistent with a uniform 
distribution of commuting X a whose eigenvalues are distributed within a ball of radius R = 2.0. 

In figure 7, we show S and the specific heat as a function of r for a = 5.0 and we observe that the 
transition occurs at r = 0.219 ± 0.005 a value consistent with r = |. We have chosen a = 5.0 to represent 
the typical behaviour for a > <3„(t = 2/9) = 4.02, so that the curve crosses the critical line r* = 2/9. 
Our theoretical prediction for the critical value of the specific heat when the transition is approached from 
the low temperature phase, eq. (|22[) predicts C v (t — 2/9, a — 5) = 1.02, which is in agreement with our 
simulations. Figure 7 also shows that 5 has a bend but there is no apparent, jump. It also appears that the 
bend occurs at a slightly larger value than r = |. The precision of our numerical results show that the bend 
occurs at r = 0.222 ± 0.002, showing small deviations when a is close (and larger) to S*(t = 2/9) = 4.02 
where the transitions take place at slightly larger values than r = § (see figure 3), however this might be 
due to finite size effects. The specific heat shows there are no strong fluctuations, so the classical theory 
should be a good approximation and in fact predicts the critical point r* = 2/9 see eq. (|12|) . The finite size 
effects in the specific heat are consistent (|31|l and account for the deviation from the limiting large N value 
of C v = 0.75. 



4.1 Eigenvalue distributions 

By measuring the eigenvalue distribution of the matrices X a and the commutators i[X a ,Xb] we can in- 
vestigate in more detail the features of the configurations at quantum level. The model has global SO (3) 
invariance, therefore each of the three matrices has the same eigenvalue distribution and also each of the 
commutators share a common distribution. We also measure the eigenvalues for the 2N x 2iV matrix 



C = O a X a 



(35) 



where a a are the three Pauli matrices. This matrix is particularly useful because it encodes information from 
the three matrices X a simultaneously and its spectrum can be easily computed for specific configurations. 
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Figure 8: Expectation value of the action < S > and specific heat C v for r = 0.23, 0.24 as a function of 
a = a/\N. These values of r lie in the matrix phase. 



For instance, when X a — aL a , where L a are IRRs of SU(2), with L a L a — 4 , we have 

2 2 1 2 JV 2 

(o- a L a ) — L a — a a L a then (a a L a + -) = — 1 



(36) 



and we find that with C —¥ U 1 Cdi ag U with U € U(2N) and noting that C is traceless and is the SU(2) 

\ _ l 

2 2 



tensor product 2 <S> N = (N + 1) © (iV — 1) we see that the eigenvalue multiplicities are JV — 1 for 



and TV + 1 for % ~ 1 therefore 



/(jv-n (-jv-i) 1 

'diag = a I ljV + 1 © ljV-1 



More generally for a reducible representation J a of ([SJ we have 
Cdiag = a I 



lrii+1 © ln,-l + U 



(37) 



(38) 



where corresponds to the one dimensional representations in J a . We see that C is sensitive to the repre- 
sentation content of the matrices. The spectrum of C will have distinct eigenvalues for the different 5(7(2) 
IRRs present in X a . 

We also measure the eigenvalues of the Dirac operator 



<7a [X a , 



(39) 



where [X a , ■] means X a acting as commutator: [X a , M] = X a M — MX a for any matrix M G MatN- 

The spectrum of D = a a [X a , ■] can be easily computed when X a = aL a - For this particular case we have 



&a[X a , ■] = a a ■ C = a 



rr2 _ r 2 _ 3 



j(j + 1 ) - 1(1 + 1) 



with j — I ± 1 with I taking the values I — 0, 1, .., JV — 1 and we obtain the spectrum 



spec{D} 



n <| -(/ + !) 





j = l+%, g(l) = 2(1 + 1) 
3 = l+h 5(0 = 2/ 

/ = 0, g(l) = 2 



(40) 



(41) 



g(l) is the degeneracy. The Dirac operator D has therefore the spectrum 

spec{D} = a{-N, -(N - 1), ... , -3, -2, 0,1,2,3, ... ,(N - 1)}. 



(42) 



This spectrum is well reproduced in our numerical simulations for the parameter range of the phase diagram 
Figure 3 corresponding to the fuzzy sphere, see Figure 10. 



4.2 Numerical results for eigenvalue distributions 

4.2.1 The t = case 

We begin with our results for the case r = 0. This situation corresponds to the model studied in [1511191 [20] 
and the system crosses the phase boundary of Figure 3 at a — 2.08 as predicted from the effective potential 
or a = 2.01 ± 0.01 as measured from simulations. 

In the low temperature phase (large 3), the dominant configurations are IRRs of SU(2), D a ~ L a . Figure 
9 shows the eigenvalue distribution for N — 24 and a — 5.00 and r = 0. The simulation indeed shows that 
the spectrum of both D3 and i\D\, D2] are discrete and equally spaced consistent with an IRR of SU(2). 

In figure 10, the eigenvalues for the matrix C = a a D a and the Dirac operator T> = a a [D a , •] are showqj 
and are consistent with the spectrum ([37p for matrix C, and (1421) for D. 
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Figure 9: Eigenvalues for D 3 and i[£>i,D 2 ]. We observe a discrete spectrum. N — 24, a — 5.00, r = 0.0. 
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Figure 10: Eigenvalues for matrix C = a a D a and Dirac operator V — a a [D a ,]. N — 24, a — 5.00, r = 0.0. The 
eigenvalues for the matrix C are nicely distributed in two peaks around C = — 4p = —12.5 and C = +^ = 11.5 
as expected. While the operator V reproducing the spectrum (|4"2l . 



In the matrix phase, we observe a continuous spectrum for the matrices X a ; we also observe that the 
eigenvalue distribution for X a has an increasing number of oscillations and becomes a smooth convex spec- 
trum for large N and as might be expected from ([23p the spectrum of X a is largely independent of a. The 



1 Note: We have defined 
this would now be T> + <j>. 



aC and 



dV. Also, the operator T> differs from the standard Dirac operator |34| on the sphere, 



eigenvalue distribution of each X a is symmetric around zero and its support is localised within the interval 
[—2, 2] and the parabola law ([25} fits excellently with R — 2.0. The distribution for the commutator is also 
symmetric around zero. 

The matrix C = a a X a shows a distribution with two maxima around ±2. The distribution, however, is 
not symmetric with the peak on the right hand slightly higher than the one on the left. There is a small 
effect of non-zero a which is larger than any effect on X a or i[X a ,Xf,]. The distribution for the operator 
D has three peaks for small N, with the central peak around zero disappearing as TV is increased. The 
distribution also shows a slight distortion with the right hand peak larger than the left. Numerical results 
are shown in figure 11 and 12 for TV = 24, 35 with a — 1.00, in figure 12 we also show the case TV = 24 with 
a = 0.60 in order to compare. 
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Figure 11: Eigenvalues for X3 and i[Xi,X2] in matrix phase. TV = 24, a — 1.00, t = 0.0. 
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Figure 12: Eigenvalues for matrix C = u a X a and Dirac operator ED = cr a [X a ,] in the matrix phase TV = 24, 
a — 0.60, 1.00, and TV = 35, a — 1.00 for r = 0.0. The eigenvalue distributions for the matrix C and operator 
D are asymmetric 

We highlight the special case of both r = a = which we refer to as the pure Yang-Mills matrix model 
and show the distributions of X a , i[X a ,Xb], in figure 13 and those of C and D in figure 14. These are all 
symmetric and consistent with our interpretation of the matrix phase as fluctuations around commuting 
matrices whose joint eigenvalue distribution is a solid ball of radius R = 2.0. 

The eigenvalue distribution for C = a a X a has two peaks at approximately ±1.9. In the large TV limit 
the Dirac operator D, has two peaks peaks located at approximately ±2.4 and support of its spectrum lies 
in the interval 1—4,41. 



4.2.2 The r ^ case 

We again focus on the eigenvalue distributions of the observables; D a , i[D a , Db], C — a a D a and the operator 
Dirac T> in the fuzzy sphere phase and X a , i[X a , Xi], C and D in the matrix phase. By a close inspection of 
the spectrum of these matrices we confirm that the numerical results are in good accord with <j> as predicted 
by equation ([9]) in the fuzzy sphere phase as we will now describe. 
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Figure 13: Eigenvalues for X$ (left) and i[Xi, X2] (right) for pure Yang-Mills matrix model for different N. The 
eigenvalue distribution of X3 is fit by the parabolic distribution (f25|) with R — 2.0 (solid line) and is consistent 
with a background of commuting matrices whose eigenvalues are uniformly distributed inside a solid 3-ball of 
radius R. 
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Figure 14: Eigenvalues for matrix C = a a X a and the Dirac operator 
matrix model. Both distributions are symmetric. 
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= a a [Xa,} for the pure Yang-Mills 



Fuzzy sphere phase. The eigenvalue distributions in the region of the parameter space in which the 
fuzzy sphere solution D a — 4>L a exists is delimited by the critical line ()12|) and defined for a > a» whenever 
r < 2/9. The effect of r is to increase the critical point a* according to <|12p . in other words, the critical 
temperature at which the transition occurs, is lowered. For a specific set of parameter N, a we can measure 
the value for <f> for a fixed value r. Knowing that in this phase the ground state is an IRR of SU(2), i.e. 
D ~ (/) diag(—s, — (s — 1), ■ • ■ , + (s — 1), +s) (with s, the spin labeling the IRR: s — ~ 1 ) we can extract 
the value of 4> measured in simulations since 4> = - eigenvalues(D)(= XM- eigenvalues(X)). For instance, in 
Figure 15(a) for N = 37, a — 5, r = 0.18 using the largest eigenvalue, the simulation gives (f> = 0.7539 while 
the prediction Q gives <f> — 0.75. 

Selecting the first positive and negative eigenvalues of D a , which should correspond to ev(Ls) = ±1, we 
obtain ev(D a )i = 0.73 and ev(D a )_i = —0.77, (see figure 15(b)). Taking the average modulus we can get 
an estimate for <j> of 4> = 0.75 while the predicted value from eq.© is (f> = 0.7537. We can now replot D a /(f> 
and the operator V>/<f> in figure 16; we can observe that the configurations are indeed around D a /<j> ~ La- 
Ova numerical results are in excellent agreement with the analytical predictions. 

Alternatively, we can extract <f> from the eigenvalue distribution of matrix C — a a D a . In the fuzzy sphere 
we have D a — <f>L a ■ We know that the eigenvalues of C are given by 



C*± 



N 



(43) 



Therefore by choosing for instance C+ the value for (j> can be extracted from the numerical results for the 
corresponding set of parameters. We obtain similar result by choosing C_. See figure 17. 



2 The model is unstable for r < 0, yet since tunneling is suppressed in the large N limit and the fuzzy sphere phase is in fact 
stable. 
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Figure 15: (a) Eigenvalues for matrix D a deep inside the fuzzy sphere region a = 5 and t = 0.18. (b) On the 
right a zoom for the first peak on the right side of zero eigenvalue of D a from which the value for <f> is measured. 
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Figure 16: Eigenvalues for matrix D a and the operator T> — a a [D ai •] deep inside the fuzzy sphere region a = 5 
and r = 0.18. The eigenvalues are divided by the measured value of <f> = 0.75 in order to show that D a /(f> are 
fluctuating around L a , irrep. SU(2) for N = 37. 



Matrix phase. In this phase as we have mentioned before, the dependence on N drops out if we con- 
centrate on X a = aD a = -¥=D a as our basic operator. We observe that the distribution for eigenvalues 
coincides, to our numerical accuracy, with that of the pure Yang-Mills model and there is no observable 
dependence on the parameters a, r and matrix size N. The data is well fit by the parabolic distribution 
(|25[) with R = 2.0. We also observe that the eigenvalue distribution of the commutator shows no observable 
deviation from that of the pure Yang-Mills model. Figure 18 shows the case for r = 0.15, a — 1 for different 
matrix size. 

In our simulations we observe that the spectrum of the matrix C = <J a X a and the operator D are not 
symmetric, but rather distorted with respect to the pure Yang-Mills case; see figure 19(a) for r = 0.15 and 
fixed value a — 1 and figure 20 (a) and (b) for C and figure 20 (c) for D. As we see from the graphs the 
asymmetry is due to the Myers term. We attribute this asymmetry to a residual dependence on the Myers 
term which drops out as -7=, noting that graphs correspond to a — 0,0.122 and a — 0.169 for a = 0,0.6 
and q = 1.0 respectively. 

The operator D has a continuous spectrum and turning on the Meyrs term again gives an observable 
asymmetry to its spectrum, see figure 19(b) for a — 1.0 and fixed value r = 0.15 and figure 20 (c) for a 
range of parameter values. It is only for r > | that the asymmetry survives the large N limit since for all 
r < I increasing JV for fixed a drives across the transition line and into the fuzzy sphere phase. 

From the distributions we can see that the transitions that takes place in the model is the one which goes 
from a non-commutative fuzzy sphere as geometrical background to a solid 3-ball with fixed radius R = 2. 



2 iV^ 
Xj = — 



1 2,2 
— a <p 



x:<4 



(44) 



Or if we rescale to the D a we get a collapse of a sphere to a point. 
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Figure 17: Here we use the eigenvalue distribution of the matrix C = a a D a in order to extract the value of <p 
since & r N _X = <P- The parameters are given by a = 5 and r — 0.18, for different values of N. We can see that 
the numerical value for <j> is in very good agreement with the theoretical prediction (6 — 0.7537 given by eq.(9). 
On the left the corresponding plot for C_ , where (f> is given by <f> 
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Figure 18: Eigenvalues for X3 and i[X\,X<2\ for different N, r = 0.15 in the matrix phase 5 = 1. Eigenvalues 
distribute uniformly inside a 3-solid ball. 



This transition occurs when a reaches the critical value a* following the critical in the phase diagram 
depicted in figure 3. The transition line is defined by the critical curve (|12|l whenever r is in the range 
< r < 2/9. The model is unstable for r < 0. However, when r < and a > «,(r), the fuzzy sphere is a 
local minimum of the model with negative action which goes to — 00 as N — > 00 and tunneling out of this 
well becomes impossible. 

In figure 21 we track the evolution of the eigenvalue distribution of X a and matrix C as a function of 
t along the line a — 5.0 starting in fuzzy sphere. We can see the fuzzy sphere disappears as we cross 



the critical value r = 2/9. We observe that r 



actually falls in the matrix phase. We find that the 



transition tracks the line r = | for different a. However, as the end point of this line is approached, i.e. 
' 4.02 we observe that in a small neighbourhood of this value the transition actually 



12 



V 107+511/17/ 

2 



occurs at r» < I see figure 3. In figure 22 we see that there appears to be some deviation from the predictions 
of ((9| where one would expect that for a = 4.02 the transition would occur at r = | ~ 0.222, however we 
find the transition at r = 0.213 ±0.002. For this transition < S > shows a clear discontinuity consistent with 
the predictions from <(9j and we expect there is still a divergent specific heat though this is more difficult to 
detect. 

Figure 23 shows the history of the expectation value of the action < S > /N 2 for iV = 24, 35, a and 
t — 0.24 from different initial configurations; cold (D a — L a ), zero and random start. When starting from 
cold configuration we observe the system goes through metastable states. All reach the same energy level 
and no decay was observed after 2 19 Monte Carlo steps or sweeps, defined as the step in which all the matrix 
elements are updated according to the Metropolis algorithm. 
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Figure 19: (a) Eigenvalues for matrix C = <r a X a . (b) D = cr Q [X a ,] in the matrix phase for r = 0.15, a = 1 



Massive Yang-Mills model Since, in the model (J23J) the limit a — removes the quadratic and Meyrs 
term we treat this the massive deformation of this case separately. 
For a = we therefore also study the model 



S[X]=NTr(--[X b ,X b 



1 - 2 + ;^ 



(45) 



We study the eigenvalue distributions for matrix X a and matrix C = a a X a for different values of the 
parameter r and matrix size N . We find that the phenomenology for this model is similar to that of r = 
described above. The eigenvalue distribution for the matrix configurations X a and small r are fit by the 
parabola (f25|) . 

We also observe that the eigenvalues for the matrix C distribute as in the pure Yang-Mills model. See 
figure 24. For values § < the matrix action is not bounded from below. 



Comment on the Massive Myers model with complex Meyrs coupling. 

let us consider the following action 

V = N Tr ( ye ofc X a [X b ,X b ] + r -X; 



For completeness 



(46) 



The model is exactly solvable (see Hoppe [24]) and is equivalent to a 2-matrix Yang-Mills model with massive 
deformation. For real a the model is not stable with no ground state. It may however be possible to localise 
configurations in a well by suppressing tunneling in the large N limit. We observe that such wells exist for 
large positive r. 

5 Conclusions 



We have performed numerical simulations of a simple two parameter 3-matrix model with energy functional 
given by ([2]). In particular we studied the phase diagram as a function of the two parameters r and a. Earlier 
studies [20] looked at the case r = and found that it exhibited an exotic phase transition at a, = (|) .It 
was argued, based on an effective potential calculation that the model should have a line of phase transitions 
dividing the (r, a) plane in two and predicting the coexistence curve. We find that for < r < | the 
coexistence curve is predicted well by the theoretical expressions. The coexistence curve asymptotes to the 
line r = § and appears to deviate slightly from the prediction of the effective potential ([9]). This special 
value of r corresponds to the value where the action is a complete square, see eq. ©• We examine the 
eigenvalue distributions of different operators. In particular we look at C = o a X a which detects the SU(2) 
representation content of the configurations X a . We find the IRR of dimension N has fluctuations around 
a lower S than any other configuration in the parameter range corresponding to the fuzzy sphere phase. 
Also, in this range the spectrum of C has the two peaks corresponding to a(j>a a L a with L a the 5(7(2) IRR 
of dimension N . In the parameter range corresponding to the matrix phase C has a continuum spectrum. 
We also study the Dirac operator D = <j a [X a , ■] and find that in the fuzzy sphere phase its spectrum is (as 
expected) a shifted, cutoff version of the commutative sphere. 

We found evidence for two distinct transition types in the emergent geometry scenario. For < r < | 
we found that, as the transition is approached from the fuzzy sphere phase with fixed r, the model has a 
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Figure 20: Effect of Myers term on the eigenvalues for matrix C = a a X a (left) and Dirac operator D = cr a [D a , ] 
(bottom) in the matrix phase. We can see that the Myers term has a small effect of the distributions, they are 
no longer symmetric as compared with the pure Yang-Mills (in black thick line). 



divergent specific heat with critical exponent a 



However for a > <5* 
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1/4 



4.02, 



and crossing the transition at fixed a, there appear to be no critical fluctuations; the transition is one 
with a continuous, but non-differentiable internal energy, and a discontinuous specific heat. In all cases the 
transition is from a fuzzy sphere to a matrix phase. We study the matrix phase in detail and find that a 
useful description of this phase is in terms of fluctuations about a background of commuting matrices whose 
eigenvalues are concentrated within a sphere of radius R = 2.0. This is consistent with the estimates of [6] 
who performed a 2-loop analysis and estimated R ~ 1.8 

We find that, though it is possible for configurations other than the irreducible fuzzy sphere to be present 
in the model, they never correspond to the true ground state of the system. Such configurations, were they 
present, would be easily detected by the matrix C. Furthermore for r = | we find that < S > when trapped 
in a fuzzy sphere configuration is larger than that for fluctuations in the matrix phase, and the fuzzy sphere 
is not a true ground state of the system. Also, for r > | we observe decays from the fuzzy sphere to the 
matrix phase, i.e. the fuzzy sphere is a meta-stable configuration for the system, with observable decay. In 
contrast for the critical line r = | we observe no decay of the fuzzy sphere to the true ground state. We 
infer that the barrier is sufficiently high in this case that the limit of N — > oo prevents tunneling out of the 
local minimum corresponding to the fuzzy sphere. 

A natural further stage in the study undertaken here is to include the effect of Fermions. However, the 
most interesting case involves either complex actions or fluctuating signs in the Fermionic sector [35] both 
of which lead to significant numerical difficulties. 

In d — 10 the pure Yang-Mills model has received significant attention in a scenario of emergent gravity 
[30] . It appears from our study, that for finite N, all such configurations will be meta stable. In the large 
N limit tunneling will be suppressed and these states may become stable, unfortunately, reliable numerical 
simulations in such a situation are more difficult. 
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Figure 21: Evolution for eigenvalues for matrix X3 for fixed 5 = 5 for different values of r. At r = 2/9 the 
fuzzy sphere disappears. On the right eigenvalues for the matrix C = a a X a . 
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Figure 22: The spectrum of < 
t = 0.213 ±0.002. 



= o a X a for varying r at <S(r = 2/9) = 4.02. We see that the transition occurs at 
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